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Abstract 

The electronic structure of a single monatomic gold wire is presented for the first 
time. It has been obtained with state-of-the-art ab-initio full-potential density- 
functional (DFT) LMTO (linearized muffin-tin orbital) calculations taking into ac- 
count relativistic effects. For stretched structures in the experimentally accessible 
range the conduction band is exactly half-filled, whereas the band structures are 
more complex for the optimized structure. By studying the total energy as a func- 
tion of unit-cell length and of a possible bond-length alternation we find that the 
system can lower its total energy by letting the bond lengths alternate leading to a 
structure containing separated dimers with bond lengths of about 2.5 A, largely in- 
dependent of the stretching. However, first for fairly large unit cells (above roughly 
7 A), is the total-energy gain upon this dimerization comparable with the energy 
costs upon stretching. We propose that this together with band-structure effects is 
the reason for the larger interatomic distances observed in recent experiments. We 



Preprint submitted to Elsevier Preprint 



1st February 2008 



find also that although spin-orbit couplings lead to significant effects on the band 
structure, the overall conclusions are not altered, and that finite Au2, AU4, and Au6 
chains possess electronic properties very similar to those of the infinite chain. 
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Very recently, single chains of suspended gold atoms have been produced be- 
tween two [110]-oriented tips in transmission electron microscope (TEM) [1] 
and in mechanically controllable break-junction (MCB) experiments [2]. In 
the TEM experiment, conductance measurements were performed simultane- 
ously with electron-microscope images of the atomic-size contact as the tips 
were separated. Before breaking, the contact consists of a monatomic nanowire 
made up of four gold atoms and shows a conductance of 2e 2 /h. The MCB ex- 
periment, although lacking direct imaging of the systems, gives also evidence 
of the formation of a monatomic wire upon stretching. 

The formation of such a monatomic chain structure for gold and its properties 
upon stretching were studied theoretically by S0rensen et al. [3] by means of 
classical molecular-dynamics simulations although for contacts between [100]- 
oriented tips. More recently, Torres et al. [4] studied the thermodynamical 
foundations of the spontaneous thinning process as well as the stability of the 
monatomic gold wire for the [110]-oriented tips using classical many-body force 
simulations as well as ab-initio local-density (LDA) and generalized-gradient 
(GGA) density-functional electronic-structure calculations. They found that 
for bond lengths above 2.8 A, the system will break into isolated Au 2 dimers. 
However, in the experimental studies, the finite monatomic wire shows an 
overwhelming stability upon stretching and is stable for bond lengths up to 
between 3 and 4 A . 

In the present work the problem of the stability of a single monatomic gold 
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wire is addressed. The electronic structure of the system is presented for the 
first time and its behaviour upon stretching used to rationalize the whole 
scenario. In particular we show that, upon stretching, the ordering of the 
energy levels leads to a situation with one exactly half-filled electronic band so 
that dimerization (i.e., bond- length alternation) becomes favoured. We show 
how, as the system is stretched, the energy gain upon dimerization increases, 
but the nearest-neighbour interatomic distance at equilibrium stays constant 
at approximately 2.5 A. Finally, the importance of relativistic effects, not 
considered in any previous study, is addressed and the results for the infinite 
chains are compared with similar ones for finite Au 2 , Au 4 , and Au 6 chains. 

We performed full-potential density-functional LMTO (linearized muffin-tin 
orbital) calculations with a local-density approximation (LDA) on isolated 
periodic infinite gold chains using the method described in Refs. [5,6]. This 
method is specifically targeted for isolated, infinite, periodic, helical polymers 
and chain compounds and have been applied successfully to a wide series 
of systems [7]. The basis set consists of two sets of s,p and d functions on 
all sites; each function is defined numerically inside non-overlapping atom- 
centered spheres and in terms of spherical Hankel functions in the interstitial 
region. The two sets differ mainly in decay constants of the latter. Scalar 
relativistic (SR) [8] corrections were included in all the calculations presented 
here; in addition the effects of spin-orbit (SO) couplings were also considered. 

For the calculations on the isolated, undimerized gold chain we assume that 
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the nuclei are lying in the (x, z) plane with the z axis parallel to the chain axis. 
We use 31 equidistant points in half-part of the first Brillouin zone in order 
to ensure the appropriate convergence of all the physically relevant quantities 
and to properly describe the metallicity of the system. We also use the level- 
broadening scheme described in [9] with an electronic temperature of 0.068 
eV (5 x 10~ 3 Ry). For the dimerized chain we use 16 points in half-part of the 
first Brillouin zone. 

The structure of the undimerized chain can be described with the single bond 
length a, whereas the dimerized chain has two alternating bond lengths a\ and 
a 2 . From these we define an average bond length 

a = ^(ai + a 2 ) (1) 

and a dimerization coordinate 

8= l -{ ai -a 2 ). (2) 

The electronic band structures of a single monatomic gold wire with an in- 
teratomic distance of a =3.5 A are shown in Fig. l.a. This bond length is 
representative of those observed experimentally. The system is metallic with 
a half-filled band of symmetry a. Below this, a broader occupied a band as 
well as (doubly degenerate) narrower n and 5 bands are found. Analyzing the 
orbitals it turns out that the a bands have contributions from both 6s and 
5d z 2 functions, whereas the n and 5 bands largely are due to d functions. 
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Including the spin-orbit (SO) couplings leads to the band structures of Fig. 
l.b. Except for the fact that the lower symmetry leads to a splitting of the 
doubly degenerate bands as well as to avoided crossings between various pairs 
of bands, the overall picture is not altered and, most important, the occurrence 
of one exactly half-filled band remains. 

Finally, Fig. l.c shows the energy levels for finite Auat chains consisting of 
N = 2, 4, and 6 atoms, respectively. In these calculations we set all bond 
lengths equal to 3.5 A, but we stress that these do not correspond to the 
optimized values (for the dimer, the optimized bond length is in fact much 
shorter, as we shall see below). Instead, they are similar to the ones observed 
experimentally. 

In Fig. l.c it can be seen that although increasing the number of atoms leads 
to some broadening of the energy regions spanned by the orbitals, most of 
the features of the infinite systems are recovered already for these fairly small 
systems. Most notably, the fact that a orbitals are those appearing closest to 
the Fermi level and that n and 5 orbitals appear at deeper energies is true also 
for the finite chains. Furthermore, the Fermi energies of the finite systems are 
very similar to those of the infinite chains. 

In Fig. 2 the total energy of the nanowire is shown as a function of the inter- 
atomic distance a. Empty squares (triangles) indicate values of the cohesive 
energy for which the relativistic contributions were included up to the SR (SO) 
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level. At the SR level, an equilibrium distance of about 2.65 A is predicted to- 
gether with a total-energy minimum of 1.35 eV/atom in good agreement with 
the results of the other ab-initio calculations [4]. The inclusion of the SO cou- 
pling leads to a reduction of the equilibrium length to 2.55 A and an increase 
of the cohesive energy by approximately 0.10 eV/atom. Such contractions due 
to relativistic effects are often observed (see, e.g., [10]). 

In both curves of Fig. 2 the energy values between 2.7 and 2.9 A are not shown. 
For these, numerical problems obscured the calculations, i.e., the highest oc- 
cupied 7r band for k = (cf. Fig. l.a) was lifted to so high energies that it 
became partly empty. This placed the Fermi level very close to a van Hove 
singularity in the density of states leading to smaller discontinuities in the 
total-energy curve. Although the effects were very tiny they were observable. 
In addition, they show how sensible the system is to external perturbations. 
Below 2.7 A the doubly degenerate tt band at k = in Fig. l.a is lifted even 
further so that this band is only partially filled and the broad a band no longer 
is exactly half-filled. 

An exactly half-filled band, as found for a above 2.9 A, favours strongly a 
(Peierls) dimerization. Therefore, in Fig. 3 the energy gain upon dimerization 
is shown as a function of the dimerization coordinate S for some selected values 
of the average bond length a that lie in the experimentally accessible region. In 
contrast to previous calculations [4], we find that the system possess a stable 
structure consisting of alternating shorter and longer bonds where the shorter 
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bonds have lengths of about 2.5 A. Furthermore, by comparing with the total- 
energy curves of Fig. 2 we see that the energy gain upon dimerization first for 
a about 3.5 A becomes comparable with the energy costs upon stretching. 

Trans polyacetylene (CH) X is a well-known example of a material that pos- 
sesses a Peierls dimerization, leading to alternating single and double bonds 
between the carbon atoms of the backbone. For this, however, both the am- 
plitude of the bond-length alternation as well as the related total-energy gain 
are much smaller than those observed here for the gold chain [6]. 

In Fig. 4 we show the band structures for two representative values of a, and for 
each case we show them for two values of the dimerization coordinate 5, i.e., a 
very small one and the optimized value in Fig. 2. The occurrence of a band gap 
at the Fermi level due to the dimerization is readily recognized for the smaller 
value of 5. For larger values of 5, at least for the larger values of a, the highest 
occupied orbital is no longer derived from the a band crossing the Fermi level 
for the undimerized structure, but from bands of ix or 5 symmetry. Moreover, 
for these larger values of a, the bands become fairly flat for the optimized 
structure which indicates that at those distances the electronic interactions 
between the dimers are only weak. For the sake of comparison we show in 
the figure also the single-particle energies for the isolated dimer with a bond 
length of 2.5 A. These are seen to lie fairly close to the band regions for 
the optimized structures supporting that this structure essentially consists of 
weakly interacting Au2 units. 
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The fact that the nature of the band gap changes from being a direct gap 
between the two a bands to becoming indirect for larger values of 5 is seen in 
Fig. 5 that shows the band gap as a function of 5 for different larger values 
of a. For smaller values of a, the above-mentioned fact that the a band no 
longer is exactly half-filled and that the n bands are partly emptied makes the 
smallest band gap vanishing up till some a-dependent threshold for 5. 

It is remarkable that for the smallest values of the dimerization coordinate 5 all 
curves lie on top of each other. Assuming that the a band crossing the Fermi 
level for the undimerized structure can be described with a single Wannier 
function per atom and that only nearest-neighbour hopping integrals need to 
be taken into account, the band gap is 

-E'gap = 2|ti — ^2 | , (3) 

where t\ and t 2 are the two hopping integrals. These may in turn be assumed 
to depend linearly on the bond lengths, 

h,2 = *o - « • (ai,2 - a). (4) 

a is an electron-phonon coupling constant. For any average bond length a one 
would expect both t and a to depend on a. However, since Eqs. (3) and (4) 
imply that 

E g&p =4a\5\, (5) 
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we obtain that the electron-phonon coupling constant is independent of a, at 
least in the range considered here. Since the tendency towards dimerization 
largely is determined by the size of the electron-phonon coupling, this result 
implies that the strength of this tendency is independent of the unit-cell length. 
Finally, the fact that the curves of Fig. 5 for larger 5 do not lie on the top of 
each other is due to the above-mentioned change of the nature of the gap. 

Experimentally, it is found that the finite chains stay stable up to bond lengths 
of about 4 A, after which the chains break. Compared with the systems we 
have studied here there is a number of differences. First, the experimentally 
studied systems consist of only about four atoms, but as Fig. 1 shows, already 
such systems possess electronic properties close to those of the infinite chain. 
Second, the experimental systems are not isolated but suspended between 
two tips, which may give further support for studying infinite chains. Third, 
and more importantly, the experimental systems are not static but subjected 
both to mechanical (stretching) forces and to electrostatic forces (due to the 
voltage between the tips). Here, we have only considered the static parts of 
the mechanical forces, which, however, due to the very different time scales 
between structural relaxations and the applied forces is justified. On the other 
hand, the applied electrostatic forces may influence the physically properties 
since they lead to an overall asymmetric potential along the chain, although 
the potential is weak (in the 10 meV range). 

In total, our results suggest the following for the experimental systems. For a 
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given structure with a certain average bond length a obtained through stretch- 
ing, the system may attempt to lower its total energy upon structural re- 
laxation. Here, our results show that the overall driving mode for this is to 
split the system into more or less strongly interacting dimers, although band- 
structure effects give that this happens first for a above around 3 A. As a 
competition to this, the system may attempt to relax towards the shorter, 
optimized, undimerized structure. Although this latter is prohibited by the 
external mechanical forces, we may still compare the two relaxation modes 
and, then, first for average bond lengths above about 3.5 A is the former 
energetically preferred. Furthermore, due to the external voltage, the dimer- 
ization mode is weakened. Therefore, we suggest that first when a is so large 
that the energy gain upon the two relaxation modes are comparable, will the 
system change structure, i.e., split into fairly well separated dimers. This of- 
fers thus an explanation for the unusually long average bond lengths that are 
observed experimentally. 

In conclusion we have shown here how it is necessary to analyze the electronic 
band structures of a single monatomic gold chain in order to gain more un- 
derstanding into the puzzling problem of its stability. In particular, we have 
shown that the dimerization is the most relevant structural relaxation that 
shall be considered for these chains. Furthermore, we have shown that rela- 
tivists effects have significant effects on the band structures although, maybe 
surprisingly, without changing the general picture. To our knowledge, rela- 
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tivistic effects have not been considered previously for such systems. Finally, 
we demonstrated that the finite chains have properties very similar to those 
of the infinite chain. 
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Figure 1. (a-b) Electronic energy bands of an infinite, monatomic gold wire and (c) 
electronic energy levels of some finite linear Au/v chains. The interatomic distance 
for all calculations is set equal to 3.5 A. In (c), +, o, and x mark ir and S orbitals, 
empty a orbitals, and filled a orbitals, respectively. In (a-b), k = and k = 1 
represent the center and edge of the first Brillouin zone, respectively, and the vertical 
dashed lines the Fermi level. 

Figure 2. The cohesive energy in eV per atom for a monatomic gold wire with 
non-alternating bond lengths. Empty squares (triangles) represent results obtained 
including relativistic contributions up to the SR (SO) level. Lines are drawn as a 
guide-to-the-eye. 

Figure 3. Energy gain upon dimerization for selected unit-cell lengths. The results 
are shown as a function of the dimerization coordinate 5, and the different curves 
are obtained for the different average Au-Au distances a shown in the figure. 

Figure 4. Band structures for infinite chains with alternating bond lengths The 
extra symbols in (b) and (d) represent the single-particle energies for an isolated 
Au dimer with a bond length of 2.5 A and using a notation as that of Fig. 1. The 
values of (a, S) in the four panels are (a) (3.0,0.1), (b) (3.0,0.4), (c) (3.6,0.1), (d) 
(3.6,1.1) in A. 

Figure 5. The smallest band gap between occupied and unoccupied orbitals as a 
function of the dimerization coordinate 5 for the same values of a and using the 
same notation as in Fig. 3. The horizontal dashed line is the value for the isolated 
Au2 dimer with a bond length of 2.5 A. 
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gure 1: L. De Maria and M. Springborg, Electronic structure. 
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